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We derive the dominant part of the 0{at 2 s ) correction to the photon energy spectrum in the 
inclusive decay B — > X s + 7. The detailed knowledge of the spectrum is important for relating the 
theoretical calculations of the B — * X s + 7 decay rate and the experimental measurements where 
a cut on the photon energy is applied. In addition, moments of the photon energy spectrum are 
used for the determination of the fo-quark mass and other fundamental parameters of heavy quark 
physics. Our calculation reduces the theoretical uncertainty associated with uncalculated higher 
orders effects and shows that, for B — > X s + 7, QCD radiative corrections to the photon energy 
spectrum are under theoretical control. 
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I. INTRODUCTION 

The process B — > X s + 7 plays an important role in 
particle physics. Within the Standard Model, it is loop- 
induced; in many extensions of the Standard Model the 
decay rate of B — » X s + 7 receives sizable contributions 
from virtual, yet undiscovered particles. Comparison of 
the experimentally measured rate for B — > X s + 7 with 
the theoretical expectations puts stringent constraints on 
various new physics scenarios. 

Thanks to BaBar, Belle and CLEO experiments, high 
quality data on B — > X s + 7 is available at present 
0, H H Due to the large irreducible background at 
small photon energies, all measurements are performed 
with a lower cut on the photon energy E 7 > E cut . Until 
recently, E cut ~ 2.0 GeV was routinely used. Last year 
Belle |l| presented precise measurements of the branching 
fraction and the first two moments of the photon energy 
spectrum with the cut as low as E cut = 1.815 GeV. This 
year, preliminary results with -E^-ut > 1.9 GeV were re- 
ported by BaBar collaboration [J ■ 

The theoretical understanding of the decay B — > X s +7 
in the Standard Model is currently quite robust. The 
total decay rate is known through next-to-lea ding or- 
der (NLO) inperturbative QCD [1 H & H El El 
13 El El EE El and the process of calculating it 
through next-to-next-to-leading order (NNLO) is under 
way 

m El El MM 

However, because the experi- 
mental measurements require a cut on E 1 , even a perfect 
knowledge of the total, fully inclusive decay rate is only 
useful if the the photon energy spectrum is understood 
well enough to relate a measurement of the branching 
fraction with E 7 > E cut to the theoretically known in- 
clusive total rate. In addition, the detailed knowledge of 
the photon energy spectrum has its own merit. Since the 
photon energy spectrum is largely insensitive to contri- 
butions from beyond the Standard Model physics, it is 



used to study strong interaction physics in heavy flavor 
decays. In particular, the value of the 6-quark mass, the 
average kinetic energy of a & quark in a £?-meson and the 
importance of the shape function for different values of 
the cut on the photon energy can be investigated. 

The theoretical description of the radiative &-decays is 
based on an effective theory approach where all particles 
heavier than the b quark are integrated out. The full 
effective Lagrangian that is used in the calculation of ra- 
diative decays of B-mesons can be found in [T2I |22| . It is 
well established [2! [24| that the dominant contribution 
to the photon energy spectrum comes from the local op- 
erator O7 ~ Sa^bF^, while for the computation of the 
total decay rate, other operators are also important. 

The photon energy spectrum in b — > X s + 7 is cur- 
rently known through 0(a s ). The effects of higher order 
QCD corrections are traditionally studied |24| in the so- 
called Brodsky-Lepage-Mackcnzic (BLM) approximation 
[25| ; the BLM corrections are both the easiest to compute 
and often provide the dominant part of the 0(a 2 s ) cor- 
rections. In flavor physics, where typical energy scales 
are relatively low and hence the strong coupling con- 
stant is large, the BLM corrections can be significant 
and their naive application could lead to inflated the- 
oretical uncertainties. A consistent application of the 
BLM corrections in i?-decays was developed in a series 
of papers [H|23| . Similar approach to the photon energy 
spectrum in B — > X s + 7 h as been initiated in and 
later elaborated upon in |29ll3fJ| where, in particular, the 
BLM corrections have been resummed to all orders in the 
strong coupling constant. Among other things, this ap- 
proach requires a careful separation of perturbative and 
non-perturbative effects and understanding the intricate 
interplay between them. 

Additional theoretical information on the shape of the 
photon energy spectrum or on the integrated spectrum 
in the presence of a lower cut on the photon energy, 
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can be obtained from the universality of the soft and 
collinear gluon radiation jsj. Such an approach can 
be used to predict large logarithms ~ ln(i? max — E cut ), 
where -E max is the maximal energy allowed for the pho- 
ton in B — > X s + 7. At present, such computations 
are performed with next-to-next-to-leading logarithmic 
accuracy (see also On general grounds it 

is clear that such an approach is applicable only for 
E C ut/E max « 1 and not for moderate or low values of 
the cut on the photon energy. 

Further improvement in the theoretical description of 
the photon energy spectrum in b — > X s +j requires O(o 2 ) 
corrections beyond the BLM approximation and it is the 
purpose of this paper to provide them. We restrict our- 
selves to the contribution generated by the operator O7 
since, as we pointed out above, this operator provides the 
leading contribution to the photon energy spectrum. 

The paper is organized as follows. In Section II we 
present the 0(a 2 s ) correction to the photon energy spec- 
trum in b — > X s + 7 and discuss its derivation. In Section 
III we study the properties of the spectrum and compare 
it with known partial results. In Section IV we consider 
the effect of our result on the evaluation of the first mo- 
ment of the photon energy spectrum and on the branch- 
ing ratio Br[_B — > X s j] in the presence of a lower cut on 
the photon energy. Finally, we present our conclusions. 



II. THE PERTURB ATIVE EVALUATION OF 
THE PHOTON SPECTRUM 

We consider the decay of an on-shell &-quark into 
a photon and any hadronic state containing a strange 
quark. We assume that this decay is mediated by the 
effective operator O7: 



o 7 = 



32tt 



(1) 



Here, is the electromagnetic field strength tensor, e = 
V 47ra, where a is the fine structure constant evaluated 
at zero momentum transfer |34[ and m is the MS 6-quark 
mass evaluated at the 6-quark pole mass m = rn(mb). 
The photon energy spectrum is parameterized in terms 
of the variable: 



we neglect all quark masses other than the mass of the b 
quark and use dimensional regularization for ultraviolet, 
infra-red and collinear singularities. Because we restrict 
ourselves to the operator O7, the photon energy spectrum 
is finite after the standard renormalization procedure is 
performed. 

At leading order in perturbation theory the kinematics 
of the decay b — > s + 7 is very simple: the photon and 
the (massless) strange quark are back-to-back with equal 
energy and momentum. Such a simple kinematics implies 
that the energy spectrum of the photon at this order is 
given by a delta-function ~ 6(1 — z). At order 0(a s ), the 
gluon radiation smears the photon energy spectrum and 
a radiative tail down to E 1 — appears. Therefore, the 
0(a 2 s ) correction to the photon energy spectrum is the 
first non-trivial QCD correction to the radiative tail. 

For the physical observables studied in this paper, we 
only need to compute the 0(a 2 s ) corrections to the pho- 
ton energy spectrum away from the maximal value of the 
photon energy E-y < E max . As a consequence, we do not 
compute the two-loop virtual corrections to b — > s + 7 
since they only contribute for 2 = 1. Therefore, to cal- 
culate the 0(a 2 ) corrections to the photon energy spec- 
trum in b — > X s + 7 for z < 1, we have to consider 
processes with up to two gluons or a quark-antiquark 
pair in addition to the strange quark and the photon in 
the final state, as well as the virtual corrections to sin- 
gle gluon emission. To perform the calculation, we use 
the techniques introduced in [35l 13(1 l37| . The idea is 
to apply the optical theorem to b — > b transition am- 
plitude in the presence of the constraint on the photon 
energy. The constraint, 8(E~ I — zi? max ), is treated as the 
on-shell condition for a "fake" particle. Then, evaluation 
of the complicated integrals over multi-particle phase- 
space, e.g. b — > s + 7 + q + q is simplified by multiloop 
integration technology [33. To solve the integration- by- 
parts identities, we use the algorithm of Ref. [3!j imple- 
mented in [iof . 

It is convenient to present the result for the photon 
energy spectrum in a normalized form, dividing the dif- 
ferential rate dY /dz by the total width for B — > X s + 7. 
This ratio is convenient because its integral over the frac- 
tion of the photon energy z is equal to one to all orders 
in the strong coupling constant a s : 



2p 7 p b 2E. 



7 



(2) 



dz = 1. 

r dz 



(3) 



where the last equality is valid in the 6-quark rest frame. 
When the photon energy spectrum is studied within per- 
turbation theory, the kinematically allowed range for the 
photon energy implies 1 < z < 1. For our calculation, 



1 Because of the Fermi motion of the heavy quark inside the B- 
meson, the photon energy spectrum extends beyond the point 
all the way up to -E max = Ms/2, where Mg is the mass 
of the B-meson. 



Note that because of Eq.©, we can restore the 5(1 — z) 
terms in the normalized photon energy spectrum without 
explicit computation of the two-loop virtual corrections. 
Our result for the photon energy spectrum reads: 

where 

F 2 = C F F (2 ' a) + CA-F (2,na) + T R n f F (2 ' cf) , (5) 
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a s is the MS coupling constant renormalized at the b- 
quark pole mass, Cf — 4/3, Ca — 3,T# = 1/2. Also, 
n/ = 4 is the number of light fermion flavors. 
The coefficient functions read: 



z + 1 



ln(l - 2) + 



7 + z-2z 2 



(6) 
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Here, Lia(x) = J dxiLi2(xi)/xi, C(3) is the Riemann 
o 

zeta-function and [ln"(l — x)/(l — x)]+ are the plus- 
distributions defined in the standard way. The con- 
stants Si take the following values S a — 1.216, S na = 
-4.795, S n{ = 49/24 + tt 2 /8 - 2C(3)/3 w 2.474. In addi- 
tion, for simplicity, we have set the renormalization and 
the factorization scales in Eq.(@J equal to the 6-quark 
pole mass mj. The total width of the radiative transi- 
tion B — ► X s + 7 is known through 0(a a ). It reads: 

r = r(°)(i + ^c F g-^)+o( a 2 )), (10) 

where = ^\V tb V*\ 2 C 2 7 (m b ) m 3 rfT 2 (m h ). 

Our result for the photon energy spectrum Eq.|@} can 
be compared to partial results that exist in the literature. 
The correction to the spectrum proportional to the num- 
ber of massless flavors n / agrees with the result in [24j . In 
addition, the terms in the photon energy spectrum that 
develop [ln n (l — z)/(l — z)]+, n = . . . 3 singularities for 
z — * 1 agree with [32]. We note that the photon energy 
spectrum computed in this paper provides an additional 
check of the constant D 2 that was extracted in |4l| from 
our previous calculation of the perturbative heavy quark 
fragmentation function |42| . 

III. ANALYSIS OF THE PERTURBATIVE 
SPECTRUM 

In this Section we discuss the properties of the pertur- 
bative photon energy spectrum Eq.Q and compare it 
with various approximations available in the literature. 
For the illustration we use a s — a s (rrib) — 0.22 for all 
plots in this Section. A choice of a s that is more appro- 
priate for observables when a cut on the photon energy 
is applied, is discussed in the next Section. 




FIG. 1: The normalized perturbative spectrum Eq.(0J in b — » 
s + 7 retaining 0{a 2 s ) (solid), BLM (dots) and 0(a s ) terms. 



An overview of the photon energy spectrum is given 
in Fig[U where the 0(a s ), the full C(a 2 ) and the BLM 
approximation to the spectrum are compared. It follows 
that when the renormalization scale is chosen equal to 
the mass of the 6-quark, the bulk of the C(a 2 ) correction 
is provided by the BLM terms and that the non-BLM 
corrections are moderate. 

As we mentioned in the Introduction, the BLM correc- 
tions require a special treatment; given their magnitude, 
the naive interpretation of Fig^ results in a large pertur- 
bative uncertainty on the photon energy spectrum. For- 
tunately, the effects of the BLM corrections have been re- 
summed to all orders in a way consistent with the Wilso- 
nian approach to the operator product expansion (OPE) 
for b — > X s + 7 j23 . Since the impact of the BLM correc- 
tions on the photon energy spectrum is well-understood, 
we can concentrate in the following analysis on the effect 
of the non-BLM corrections. This is the new contribution 
to the photon energy spectrum derived in this paper. 
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FIG. 2: The non-BLM 0{o? a ) correction to the photon energy 
spectrum relative to the 0(a s ) correction. 
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FIG. 3: The non-BLM 0(o? s ) correction (solid), its z -> 1 
approximation (dots-dashes) , the z —* 1 approximation minus 
its value for z = (dots) and the z — * 1 asymptotics times z 3 
(dashes). 

The magnitude of the non-BLM part of the 0(a 2 s ) cor- 
rection is shown in Fig|2 where it is plotted relative to 
the photon energy spectrum through 0(a s ). The relative 
size of the non-BLM 0(a 2 s ) correction is within ±10%, 
approximately. 

It is interesting to check if the terms [m"(l — z)/(\ — 
z)] + , n = ... 3 in Eq.Q that are singular in z — > 1 limit, 
provide a good approximation to non-BLM corrections; 
Figs l3l4l illustrate this. We see that the singular terms do 
not furnish a good approximation to the exact result for 
moderately large values of z; for example, for z = 0.8, the 
z — > 1 approximation overestimates the exact result by a 
factor of three. This comparison suggests that based on 
the previously known large z behavior of the non-BLM 
piece |32| . one could not have computed its contribution 
to observables, unless very high value of the cut on the 
photon energy is in place. 

Although the z — > 1 approximation alone is insufficient 
to adequately approximate the non-BLM correction, it is 
possible to improve the quality of the approximation by 



0.4 0.6 0.8 

z 

FIG. 4: Same as FigEl Each curve is normalized to the non- 
BLM 0(a 2 s ) correction. 

combining the z — > 1 approximation with the z — > 
behavior of the photon energy spectrum. It is easy to see 
that for O7 operator the photon energy spectrum should 
vanish as z 3 , for z — > 0. Hence, we can modify the z — > 1 
asymptotics by multiplying it by z 3 . This improves the 
quality of the approximation; for z — 0.8, the modified 
z — ► 1 approximation overestimates the exact result by 
a factor 1.3. However, the best description of the exact 
result is obtained by subtracting from z — > 1 asymptotics 
its value at z — 0. While this approximation is hard to 
justify theoretically, it gives quite an accurate description 
of the spectrum for moderately large values of z. 



IV. APPLICATIONS 

In this Section, we study the numerical implications 
of the 0(a 2 ) correction to the photon energy spectrum 
derived in this paper. As noted in the Introduction, high- 
quality experimental results on B — > A s +7 become avail- 
able from the B-factories. In particular the branching 
fraction of B — > X s + 7 and the average energy of the 
photon in the presence of a cut E y > E cnt are measured 
precisely. 

The applicability of the pure perturbative description 
of the photon energy spectrum depends upon the exact 
value of E cut which, in turn, defines the degree of "in- 
clusiveness" of the process. If 2£ cu t is high, the usual 
OPE breaks down and the resummation of leading twist 
effects in the OPE leads to the appearance of the non- 
perturbative shape function. By decreasing E cut it be- 
comes possible to describe the non-perturbative compo- 
nent of the photon energy spectrum by means of the lo- 
cal OPE. It is generally believed that for E cut as low as 
1.8 GeV, we are already in the OPE regime. Fortunately, 
Belle published results on the branching fraction and 
the first two moments for precisely this value of the cut 
on the photon energy E 1 > 1.815 GeV. 

To proceed further, we have to specify the numeri- 
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cal value of the fr-quark mass. Although our result for 
the spectrum is derived in the pole scheme, it is clear 
on general grounds that the pole mass can not appear 
in short-distance quantities; instead, a consistent appli- 
cation of OPE leads to short-distance, low-scale quark 
masses that enter the pcrturbative calculations. In prin- 
ciple, the transformation from one scheme to another 
should be done in a self-consistent way, including the 
proper treatment of non-perturbative effects. However, 
as we will see, the contribution of the non-BLM 0(a 2 s ) 
correction to observables is relatively small and, there- 
fore, a simple estimate suffices. For this reason, in- 
stead of the pole 6-quark mass we use in our formulas 
nib = 4.6 GcV which roughly approximates the kinetic 
&-quark mass normalized at /i = 1 GeV • For the two 
applications that we consider below, the exact value of 
the b quark mass mostly affects the z integration region, 
through z > z cut — 2E cu t/mi,. In addition, we note that 
the scale at which the strong coupling constant has to be 
evaluated in the numerical estimates has to be smaller 
than the mass of the &-quark. This follows from the fact 
that, for a given value of E cui , the maximal value of the 
invariant mass of the hadronic system in b — > X s + 7 
is y/m b (m b - 2E cut ) ~ 2 GeV for E cut ~ 2 GeV. This 
implies that for typical values of the cut on the photon 
energy the appropriate value for the scale at which a s 
should be evaluated is fj, ~ 1.5 — 2 GeV. Having made 
these preliminary remarks, we are ready to present our 
numerical estimates. 

First, we consider the fraction of events that contain 
a photon with energy above the cut -E C ut- The second 
order QCD corrections shift this fraction by: 

SR(E cut ) = - A R (z cut ), 

A R (z cut ) = C F J dzF (2 \z), (11) 


where we have used Eq.© to express the integral over 
z > z cut through the lower part of the photon energy 
spectrum. As we explained previously, the BLM cor- 
rections are taken into account in the estimates for the 
fraction of events that exist in the literature. There- 
fore, we disregard the BLM corrections and present the 
contribution of the non-BLM corrections to Ea. HllI) in 
FigEl With our choice of the 6-quark mass and for 
E cut = 1.8 GeV, we find z cut = 0.78 and A R = 0.50 
which, depending on the value of a s — 0.3-70.35, trans- 
lates into: 



non-BLM 



-(5.5 ±0.5) x 10" 



The predictions for the fraction of events with the pho- 
ton energy larger than 1.8 GeV that are available in the 
literature are: R = 0.952±£j£| ||, R = 0.9581^029 0; 
R = 0.95 ± 0.01 H3 and R = 0.89±g;g| 0, where in 
the last result only perturbative uncertainty is displayed. 
We see that the second order non-BLM QCD corrections 



FIG. 5: The function Ar(z) Eq. lllH for an experimentally 
relevant range of z cut . Only the non-BLM corrections are 
used in the calculation. 



change the fraction of events with E 7 > E cnt by less than 
a percent and are about half of the uncertainty assigned 
to R(E cut ) in H3- 

The next observable we consider is the first moment 
of the photon energy spectrum (E 7 ) for E 1 > E cut . It 
is straightforward to derive the correction to the average 
photon energy when the second order QCD corrections 
are applied. The result reads: 



6(E, ( ) = -^(^) 2 A E (z cut ), 



A E (z cut ) = C F J dz{l-z)F^{z) 

■Scut 

1 z cut 

+C 2 F J dz(l - z)F^(z) J dyF^(y). (12) 

Zcut 

The last term in Ea. l|12|l comes from the interference of 
0(a s ) corrections to the numerator and denominator in 
the expression for (E^). We disregard the BLM correc- 
tions and compute as a function of z cnt - The result is 
shown in Fig|S] For E cut = 1.8 GeV, we get z cut = 0.78 
and obtain A E = —0.68. This leads to: 



6(E 1 



1— BLM 



(1.8 ±0.3) x 10~ 2 GeV, 



(13) 



depending on the value of the strong coupling constant. 

The current most accurate measurement of the average 
energy in B — > X s + 7 comes from Belle collaboration [l| 
and reads (£ 7 )_e t >i. 8 gov = (2.292±0.026±0.034) GeV. 
The theoretical estimates for this observable that are 
available in the literature are: (E 7 ) = 2.27t '; l y GeV 
|44|. where only the uncertainty associated with uncal- 
culated higher order corrections is displayed, and (Ej) = 
2.312 GeV |2j| (no explicit estimate of the uncertainty 
has been given in that reference). Our result Eq. i|13|) 
shows that the central values in these predictions should 
be shifted approximately by one percent. In addition, as 
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FIG. 6: The function Ae(z) Ea, l|12|l for an experimentally 
relevant range of z cu t- Only the non-BLM corrections are 
used in the calculation. 



the result of our calculation, the error bars for (£L) as- 
sociated with higher order QCD corrections given in 0] 
should be reduced to at most a percent. 

We note that although a percent shift in the central 
value for {Ej) is small, it can not be considered negligible 
given the current precision achieved in i?-physics. In par- 
ticular, the average photon energy plays a central role in 
constraining the mass of the 6-quark from B-decays. The 
correction that we have computed in Ea. (|13[) results in 
approximately —30 MeV shift in the short-distance low- 
scale mass. The value of this shift can not be neglected 
because it roughly equals the uncertainty in the value of 
the 6-quark mass that was extracted from recent fits to 
semi-leptonic and radiative i?-decays |45j . 



V. CONCLUSIONS 

The knowledge of the photon energy spectrum in B — > 
X s + 7 is required to relate the experimental measure- 



ments of the branching fraction Br[B — > X s + 7] with a 
lower cut on the photon energy and the fully inclusive 
theoretical calculations. This observable is of particu- 
lar interest because it is sensitive to physics beyond the 
Standard Model and represents the cleanest observable 
for theoretical computations. In addition, the photon en- 
ergy spectrum is used to test our understanding of strong 
interactions and to extract fundamental parameters of 
heavy flavor physics. 

In this paper the analytic computation of the 0(a 2 s ) 
corrections to the photon energy spectrum is presented. 
We have restricted ourselves to the contribution of the 
operator O7. Partial analytical results on the 0(a 2 ) cor- 
rections to the photon energy spectrum that have been 
published previously are confirmed. On the other hand, 
we have shown that the terms that are singular in z — > 1 
limit do not furnish a good numerical approximation to 
0{a 2 s ) corrections for z < 0.95. 

To evaluate the numerical impact of the 0(a 2 ) cor- 
rections, we separate the BLM and the non-BLM pieces. 
While the BLM corrections have been computed and ana- 
lyzed previously, the non-BLM component at O(ot 2 ) rep- 
resents the new result of this paper. The impact of this 
new result on the fraction of events with photons having 
energies larger than E cut and on the average energy of 
the photons in B — > X s + 7 is in the 1% range. This new 
contribution is comparable in size to both the accuracy 
of the experimental measurements and the uncertainty 
assigned to these observables in previous evaluations. 
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